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rwere  represented  as  elastic  and  temperature  dependent. 
Numerical  simulations  of  a  series  of  isothermal  creep  tests 
were  performed  to  determine  th  effective  creep  compliance 
parameters  of  the  composite  constitutive  equations.  The 
macromechanical  response  of  a  composite  structural  element, 
as  predicted  by  the  effective  constitutive  equations  and 
their  derived  parameters,  was  then  verified  using  results 
computed  using  a  micromechanical  model  which  explicitly 
included  the  fiber  and  matrix  as  discrete  phases.  To 
determine  the  potential  existence  and  form  of  a  composite 
complex  modulus,  the  response  of  unidirectional  composite 
structural  elements  to  simultaneous  sinusoidal  temperature 
and  mechanical  loads  was  investigated.  Finally,  solutions 
of  free  vibration  and  transient  dynamic  analyses  of  some 
simple  composite  structures  were  performed  to  examine  the 
effects  of  the  thermoviscoelastic  behavior  on  the  damped 
response  of  some  simple  composite  structures. 
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1 .  INTRODUCTION 


This  report  summarizes  the  work  associated  with  the  evaluation  of 
the  time  and  temperature -dependent  response  of  unidirectional  fiber 
reinforced  composite  structures  and  structural  elements.  An  examination  of 
this  response  is  of  fundamental  importance  in  understanding  the  behavior  of 
fiber  composite  structures  subjected  to  an  environment  which  includes  severe 
temperature  variations.  The  proposed  large  orbiting  space  structures,  e.g. 
solar  electric  power  stations,  communication  satellites  with  large  antennae, 
etc.,  are  examples  of  structures  which  will  undergo  severe  temperature 
cycling.  Thus,  the  utilization  of  graphite/polymer  composites  in  these 
designs  requires  that  the  thermo-mechanical  response  of  structural  elements 
composed  of  such  materials  be  investigated. 

The  carbon  fiber/polymeric  matrix  composites  are  known  to  exhibit 
time  dependent,  or  viscoelastic  behavior.  In  addition,  the  time  dependent 
characteristics  are  modified  by  the  temperature  environment  experienced  by 
the  material,  so  that  a  thermoviscoelastic  representation  must  be  used.  In 
particular,  a  thermorheologically  complex  material  (TCM)  model  of  the 
matrix,  and  thus  the  unidirectional  fiber  composite,  has  been  employed  here. 

The  objectives  of  the  current  efforts  in  the  thermoviscoelastic 
analysis  of  fiber  reinforced  composite  structures  were  threefold.  Initially, 
attention  was  focused  at  developing  the  form  and  parameters  of  the 
thermorheologically  complex  constitutive  relations  for  unidirectional  fiber 
composites  of  the  polymeric  matrix  type.  Following  the  successful 
characterization  of  the  effective  constitutive  relations,  the  response  of 
unidirectional  TCM  composite  structural  elements  to  simultaneous  sinusoidal 
temperature  and  mechanical  loads  was  determined,  in  order  to  investigate  the 
existence  of  complex  moduli/compliances  in  this  situation.  Finally, 
temperature -dependent  eigensolutions ,  and  dynamic  analyses  of  transient  load 
and  temperature  conditions,  on  some  simple  composite  structures  were 
performed,  with  the  purpose  of  investigating  the  extent  and  sources  of 
damped  response  in  the  TCM  composite  structures. 


THERMOVISCOELASTIC  CONSTITUTIVE  RELATIONS 


2.1.  DEVELOPMENT  OF  EFFECTIVE  THERMOVISCOELASTIC  EQUATIONS 

The  effective  constitutive  equations  for  unidirectional  composites 
are  dependent  on  the  mechanical  behavior  of  the  individual  composite 
constituents.  In  this  work,  the  fibers  were  represented  as  transversely 
isotropic  and  linearly  elastic,  temperature  dependent  elements.  The  matrix 
material  surrounding  the  fibers  was  taken  as  isotropic.  The  deviatoric  or 
shear  components  of  the  matrix  response  were  assumed  to  be  linear 
viscoelastic  and  thermorheologically  complex.  The  dilational  components 
have  been  considered  to  be  elastic  and  temperature  dependent.  Note  that 
this  characterization  of  matrix  behavior  as  elastic  in  dilatation  and 
viscoelastic  in  shear  is  consistent  with  the  typically  assumed  response  in 
three-dimensional  viscoelastic  stress  analysis  [1]  . 

The  development  of  the  effective  constitutive  equations  was  performed 
by  assuming  that  the  form  of  the  composite  stress-strain  equations  in  all 
deformation  modes  (dilational  and  deviatoric)  was  identical  to  the  form  of 
the  matrix  constitutive  equation  in  shear.  For  typical  fiber  volume 
fractions,  this  is  a  reasonable  assumption.  The  most  simple  form  of  a 
thermorheologically  complex  viscoelastic  stress-strain  equation  was  used; 
this  is  one  which  includes  a  horizontal  shift,  representing  time -temperature 
equivalence,  and  a  single  vertical  shift  for  modeling  the  temperature 
dependence  of  the  initial  modulus.  This  is  a  slight  variation  of  the 
thermorheologically  complex  equation  proposed  by  Schapery  [2] . 

The  form  of  the  thermorheologically  complex  viscoelastic  material 
model  representing  the  matrix  material  in  shear  is  given  by  the  symbolic 
equation: 

2e.j  “  g  ®  s„  (2.1.1) 

which  represents  the  more  complete  form 

ds . . ( t ' ) 

2e  (t)  -  g(0)V (4>)  s  (t)  +  Ag(f)  s  (0)  +  f*  Ag(f-f’)  •  df  (2.1.2) 

J  J  J  n 
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where 


t* 


eij’sij  "  the  deviatoric  strain  and  stress  respectively, 


g(0)  -  the  instantaneous  shear  creep  compliance, 


g(£)  -  the  time  dependent  shear  creep  compliance, 


V(<£)  -  a  vertical  shift  factor,  and 


5 
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m 

■v 
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<i>(t )  -  temperature. 

Also,  the  parameter  £  represents  the  reduced  time  and  is  given  by  the 
equal  ion 


£(t)  -  J*^h[^(u)  ]du 


(2.1.3) 


where 


h(^)  -  a  horizontal  shift  factor. 


■rj; 


Note  that  the  temperature  dependence  of  the  stress-strain  equation 
is  represented  by  the  horizontal  and  vertical  shift  factors.  The  reduced 
time  allows  the  viscoelastic  response  at  any  temperature  to  be  obtained  from 
a  master  response  at  reference  temperature  through  a  horizontal  shift  along 
the  time  axis.  The  role  of  the  vertical  shift  factor  is  to  model  the 
temperature  dependence  of  the  initial  compliance.  Note  that  in  isothermal, 
temperature  independent  viscoelasticity,  both  the  horizontal  and  vertical 
shift  factors  are  unity.  Further  background  into  the  development  of  the  TCM 
constitutive  equations  may  be  found  in  Schapery  [2]. 

With  respect  to  the  thermoviscoelastic  unidirectional  fiber 
composite,  since  the  matrix  material  is  transversely  isotropic,  and  because 
of  the  random  nature  of  fiber  placement,  the  composite  is  statistically 
transversely  isotropic.  Its  effective  thermoviscoelastic  constitutive 


relations  may  be  written  using  the  simplified  notation  of  equation  (2.1.1) 


as 

"ll  “  S11  ®  *11  +  S12  ®  *22  +  §12  ®  *33  +  aL  ®  * 

*22  -  S*2  0  au  +  S2*  0  a22  +  S2*  ®  *33  +  ©  * 

e33  —  S^2  0  <7^  +  S23  ©  a22  +  ^22  ©  CT33  +  aT  ©  ^  (2.1.4) 

2£12  “  S44  ©  *12 

2*13  “  S44  ©  *13 

2«23  _  S55  ©  *23 

Where 

S.j  -  the  effective  creep  compliances,  and 

"fc  lit 

q^,  a 3,  -  the  effective  compliances  for  the  longitudinal  and  transverse 

coefficients  of  thermal  expansion. 

Note  that  the  bars  over  the  strains  and  stresses  indicate  average  strains 
and  stresses  in  the  composite. 

The  complete  set  of  effective  creep  compliance  parameters  include 


ft 


■f. 


-  Che  effective  instantaneous  creep  compliances, 

-  the  effective  time  dependent  creep  compliances, 

-  the  effective  horizontal  shifts,  and 

-  the  effective  vertical  shifts. 

Thus,  unlike  the  case  with  elastic  materials  in  which  effective  moduli  or 
compliances  are  determined  from  the  constant  response  to  a  constant 
mechanical  or  thermal  input,  a  series  of  isothermal  creep  tests  must  be 
performed  in  each  mode  of  deformation,  to  determine  the  time  and  temperature 
dependence  for  the  entire  set  of  effective  creep  compliances  of  the 
composite  constitutive  equations. 

Consider  first  the  case  of  longitudinal  shear.  The  effective 
stress- strain  equation  is  given  by 

2i12(t)  -  S4*(0)  V4*(*)a12(t)  +  AS4*(£Vl2(0)  (2.1.5) 


S.*(0) 

IJ 


AS  *(t) 
ij 


h .*(*) 

ij 


Vij(^) 


+  sl  as4J«-D 


da12(t' ) 
dt' 


dt' 


where 


£  -  l(u)  -  JU  h44 [4>(t)  ]  dr 
-  00 


(2.1.6) 


Suppose  we  simulate  an  isothermal  creep  test  in  which 


5 


v.v 


(2.1.7) 


<t>{ t)  - 


where  H(t)  is  the  heavy-side  (unit)  step  function  defined  such  that 


H(t) 


{? 


t  <  0 
t  >  0 


(2.1.8) 


The  stress-strain  equation  (2.1.5)  gives 


2^12<t)  -  S4*(0)V4*(^)ao  +  AS4*(h4*(^0).t)ac 


(2.1.9) 


The  shift  functions  h44(^>)  and  V 44(^)  are  defined  so  that 


hJ> R>  -  R>  -  1 


(2.1.10) 


where  is  the  reference  temperature.  In  this  way,  by  performing  the 


isothermal  creep  test  at  <f>  ~  the  creep  function  S44  (t)  is  determined: 


+  AS44^c^ 


^ 6 12 ( t ’ ^r) 


(2.1.11) 


At  t-0 ,  AS44  (0)  -  0  by  definition,  so  that 


.*■-  »'•  . **  *"  *  *'•  •"“>**  »"*  «**  -  •  ♦  *  »  %*■  * *  *  “J* 


6 


(2.1.12) 


(0)  - 


2^12<0’*R> 

o  o 


and  for  t  h  0  we  find 


*s44(t)  - 


2h2( t’V  -  2;12(q.»r) 


(2.1 


Now,  the  isothermal  creep  test  at  <f>  -  <j>Q  can  be  used  to  determine  both 
h^^(<^)  and  V^(<£)  as  follows.  At  t  -  0  we  get 


<*o>  - 


2‘l2<°’»o> 

S4t(0)*ffo 


(2.1 


■jJ* 

which  provides  the  vertical  shift  V44(^Q) •  This  leaves  the  horizontal 
h^(<£o)  as  the  only  other  unknown.  At  some  other  time  during  the  creep 

A 

test,  e.g.  t  —  t,  we  get 


AS44(h44(^o)*t)  “ 


2e12<t^o>-2i12(°^o> 


AS44(* 


44 ; 


(2. 


from  which 


h44«V  -  W*  (2-X 

Thus,  by  performing  a  series  of  isothermal  creep  tests  at  a  range  of 

★  * 

temperatures,  both  h^(<£)  and  V^( <f> )  can  then  be  determined. 


.13) 


•  14) 


shift 


1.15) 


.16) 


7 


Consider  now  the  case  of  applied  longitudinal  stress  coupled  with  a 
uniform  temperature  distribution,  with  all  other  stresses  zero.  Again 
simulating  an  isothermal  creep  test,  we  apply  the  loads 


^U<t)  -  H(t)oL 


(2.1.17) 


*<t)  -  H(tU 


in  which  case  from  the  stress-strain  relations  (2.1.4)  we  have  the  two 


equations 


<ll(tlV  “  Sll(O)Vll(0o)£7L+ASll(°aL+Q[L(,:^o)*^o 


e'22(t^o)  "  S12(0)V12Wo)aL+AS12(€>‘,L+aL(t^o)^o 


(2.1.18) 


(2.1.19) 


Now  define 


«ii(t^o)  -  «ii(t^o)’Qi(t^o)*^o  Ci-1,2)  (2.1.20) 


'ft  'Jc 

where  and  ,  and  the  subscript  i  is  not  summed.  Then  from 

(2.1.18),  at  <f>  -  <f>  we  get 


Sll<°)  “  a 


gll<°'V 


(2.1.21) 


A 


ASn(t)  - 


<n<t,V-en(0’V 


(2.1.22) 


Then  at  <f>  -  <b  , 
o 


<*o> 


•n<°'+a> 

S1*(°)-L 


(2.1.23) 


A 

and  at  some  other  time  t  -  t, 


^ll^lW^  ” 


cll(t'^o)‘ell(0’^n) 


-  ASn(en) 


(2.1.24) 


from  which 


<*o> 


From  (2.1.19) 


S^O) 


622(0,<^r) 

aL 


(2.1.25) 


(2.1.26) 


and 


ASX*  (t) 


e22(t,^R)-t22(0'^R) 


(2.1.27) 


At  4>  -  4>q, 


9 


i 


and  at  some  other  time  t  -  t, 


iS12<h12  <♦„>•*>  -  - ^ - -iS12<fl2> 


(2.1.29) 


from  which 


h12<V  “  ^12//t 


(2.1.30) 


Now  consider  the  case  of  an  applied  transverse  isotropic  stress 


aT  "  ff22  "  *33 


(2.1.31) 


coupled  with  a  uniform  temperature  distribution,  with  all  other  stresses 
zero.  Again  simulating  an  isothermal  creep  test,  we  apply  the  loads 


(*22^  ~  *33^)  ~  H(t^T 


(2.1.32) 


*(t)  -  H(t)4 


Using  the  first  of  the  stress-strain  relations  (2.1.4),  the  creep 

compliance  parameters  could  again  be  generated,  but  this  will  not  be 
repeated  here.  From  either  the  second  or  the  third  of  equations  (2.1.4),  v 


sssswHr 


a 


‘t(C**o>  “  S22(0)V22(VVS23(0)V23(*o)oT 


+^S22(^22)aT+AS23(|23)crT 


kt  <j>  -  we  have 


★  ★ 


At  t  -  0 ,  AS22  (0)  -  AS23  (0)  -  0  so  that 


S22(0)+S23(0)  "  £T(0l*R)/aT 


and  for  t  +  0, 


AS2*(t)+AS2*(t)  - 


fT(t’<V~eT(0,'*R) 


For  <f>  -  <t>  and  t  »  0, 
o  ’ 


S22<0)V22^o)+S23(0)V23^o>-£T^°-V/ffT 


and  at  <f>  =  4>q  ,  t  -  t  we  get 


AS22(h2*<4o)*t)+AS23(h23(*o).t)  - 


eT(t'<^o)'£T(0’^o) 


(2.1.33) 


S2>(0)+S23(0)+AS22(t)+AS23(t)  “  £T<t>V/ffT  (2.1.34) 


(2.1.35) 


(2.1.36) 


(2.1.37) 


(2.1.38) 


Now,  from  transverse  shear  isothermal  creep  tests  we  get  the  equations 


vr/’tw,!! 


Fr#r«v#rfr#r#r«nvi 


S55(0)  “  S22(0)  ‘  S23(0)  “  2«23<0,*R)/a23 


AS  \  AS  AS  *m  2t23(t’V  '  2£23(0,i^R) 

As55(t)  -  AS22(t)  -  AS23(t)  -  g ^ 


(2.1.39) 


(2.1.40) 


For  <f>  -  <l>  and  t  -  0, 
o 


S22  <°>V22(V-S23(0)V23<V  “  2i23(0^o)ff23 


(2.1.41) 


and  for  t  -  t, 


AS22(h22^o)‘t)-AS23(h23^o)*t) 


2t23^t’^o^'2e23^0,^R^ 


(2.1.42) 


From  (2.1.35)  +  (2.1.39), 


*  5T(0,*  )  2i?3(0-V 

S22(0)  -  1/2  ■  T  a  R  +  — “ - ^ 

T  a23 


(2.1.43) 


From  (2.1.35)  -  (2.1.39), 


c(O,0)  2e„.(0,^  ) 

S  *(0)  -  1/2  •  -X-— R  -  - L- 

23  aT  o23 


(2.1.44) 


From  (2.1.36)  +  (2.1.40), 


AS„(t)-l/2 


2e23(t,^R)-2e23(0,<^R) 


(2.1.45) 
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From  (2.1.36)  -  (2.1.40), 
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AS23(t)-l/2 


fT^t’^R^'gT^0,^R^  2e23('t’^R)'2e22,<'°’^R) 
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From  (2.1.37)  +  (2.1.41), 


V22<V' 


2S22(V 


flgjg)  ^  2£23(°^R> 


'23 


and  (2.1.37)  -  (2.1.41)  gives, 


V23(V 


2S„*(*0) 


'23' 


2e23(0'^ 
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Then  (2.1.3 8)  +  (2.1.42)  results  in 


S22(h22(V^“1/2 


er(t,*o)-<T(0,*o)+  2e23(t^o)-2i23(0,^) 


23 


”AS22(^22) 


giving 


h22^o')",'’22^t 
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(2.1.46) 


(2.1.47) 


(2.1.48) 


(2.1.49) 


(2.1.50) 
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From  (2.1.38)  -  (2.1.42)  we  get 


AS23(h23(*o)-t)-l/2 


£T(t^o)'£T(0^o)  2£23(t^o)'2i23(°^o) 
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(2.1.51) 


-  as23(£23) 

from  which 


h23(V-*23^  (2.1.52) 

Now  consider  the  case  where  thermal  loads  alone  act  on  the  composite 
cylinder.  In  the  case  of  a  temperature  step  function  rise  A<£  relative  to 
the  reference  temperature  <f> the  effective  longitudinal  secant  expansion 

coefficient  is  defined  by  while  the  transverse  expansion  coefficient 

qt  is  defined  by  the  radial  surface  displacement  per  unit  radius.  Thus 


(2.1.53) 


A<f>)  -  eT(t)/A<t> 

Schapery  [2]  proposes  that  for  the  thermorheologically  simple 
material  (TSM) ,  in  the  absence  of  mechanical  loads, 


ft  '  dA^ ( r ) 

£’!“J  aL^L'^L)  dr  dr 
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(2.1.54) 


(2.1.55) 


rt  * , .  .  . dAi(r)  , 

-j  qt(4t-^)  T1— 1 ’dr 


£22"£33  J  “T^T  ‘‘T 
-  00 


dr 


where 


hL[^<T)]dr 


(2.1.56) 


and 


V  J*  V*(r)]dr 


(2.1.57) 


For  the  most  simple  of  the  TCM  forms,  it  is  proposed  that  for  a  general 
temperature  history,  again  in  the  absence  of  applied  stress  for  convenience, 


«n(t)  -  a*(0)V*(^)A*(t)  +  Aa*(eL)A<£(0) 


(2.1.58) 


o+ 


i 


k&2 


"22(t)  "  *33(t)  “  +  Aa*(^x)A^(0) 


(2.1.59) 


,dM, 
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For  a  temperature  step  function  rise 


A0(t)  -  H(t)  A <j> 


(2.1.60) 
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we  obtain  the  strains 


"ll(t)  ”  {aL(0)VL(^)+AaL(^L)}A4i  (2.1.61) 

iT(t)  -  ;22(t)  -  i33(t)  -  <[a*(O)V*(0)+Aa*(^T)j-A^  (2.1.62) 

At  t  -  0,  Aa^  (0)  -  Aa*  (0)  -  0  by  definition,  so 

*  £ i i (°)  *  *  «x(0) 

4(0)VLW  -  -t(0)VtW  -  -Jj-  (2.1.63) 

We  now  define  an  arbitrary  temperature  <f>  —  <j>  as  the  effective  thermal 

expansion  coefficients'  reference  temperature.  Since  A$  -  0  at  the 
reference  temperature  defined  for  the  composite  effective  mechanical  creep 
compliances,  some  other  temperature  would  be  more  convenient  to  use  for  the 
composite  effective  thermal  expansion  coefficients.  Since 

Ra>  -  V*R*>  “  -  h>R*>  -  1  (2-1-64) 


there  results 


*  « 

aL(°)  - 


£  11  (°  i  ^Rq,) 

A  4> 


f 1 1 ( 0 ■ ^rq ) 

^Ra'^R 


(2.1.65) 


and  similarly 


<4(°)  - 


£T(0'^R al 

^Ra^R 


(2.1.66) 


Then  at  temperature  <£q  +  the  vertical  shifts  are  computed: 


(2.1.67) 
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a*(0)A<£ 


*  £t(0-*o> 

W  "  - 


aT(0)A^ 


where  A<f>  -  4>q  - 


Again  considering 


£ll(t) 


-{a 


*(0)  +  Aa*(t)  ^  A <)> 


>} 


«T(t)  -  ^a*(0)  +  Aa*( t)  f-  A<f> 


>} 


from  which  the  effective  delta  compliances  are  defined: 


A  *,*.v  gll(t,'^Ra)  * 
AaL(t)  “  - 


A  <f> 


aL(°) 


Aa*(t) 


£T(t’V) 


A <f>  ‘  aT(0^ 


Then,  at  <f>  -  4>  and  t 
o 


t  f*  0, 


*  *  -  ;n<t,*  )-a*(o)v*c*  )  * 

toL<hL(^)-t)  -  ^ - °  J-  ^  -  toL«L> 


Aa*(h*(^o)*t) 


«T(t,^0)  -  q*(0)V*(^o) 


A  <f> 


Aa^(^T) 
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(2.1.68) 


(2.1.69) 


(2.1.70) 


from  which 


i: 

V 
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b 


£ 


hL<V  -  V‘ 


(2.1.71) 


and 


W  “  VC 


(2.1.72) 


2.2.  LITERATURE  SURVEY  OF  THERMOVISCOELASTIC  MATERIAL  PROPERTIES 


As  part  of  this  work,  an  extensive  literature  search  was  conducted 
in  an  effort  to  determine  realistic  material  properties  of 
thermoviscoelastic  solids,  in  particular  those  typically  employed  as  epoxy 
matrices  in  composites.  Several  papers,  e.g.  [3-8],  were  found  which 
provided  some  or  all  of  the  material  parameters  for  a  thermorheologically 
simple  characterization  of  typical  matrix  material.  However,  very  little 
data  exists  for  thermorheologically  complex  viscoelastic  epoxies. 

A  number  of  papers  [9-13]  and  a  private  communication  [14]  were 
obtained  which  provided  useful  data  on  thermorheologically  complex  epoxy 
matrices.  All  of  the  papers  [9-13]  used  a  characterization  scheme  which 
fitted  the  data  to  a  semi-empirical  model  first  developed  by  Schapery  [15] 
for  nonlinear  isothermal  viscoelasticity.  The  data  in  all  of  these  papers, 
first  reported  in  [9]  and  [10],  was  modified  for  use  in  the  work  reported 
here,  as  described  in  a  previous  report  [16]. 

The  matrix  thermoviscoelastic  data  used  in  this  study,  including  the 
time -dependent  matrix  shear  creep  compliance,  the  matrix  shear  horizontal 
shift  and  vertical  shift,  are  included  in  Figures  2.1,  2.2  and  2.3, 
respectively.  The  fiber  elastic  properties  and  the  matrix  properties  at 
reference  temperature  are  listed  in  Table  1. 


2.3.  EVALUATION  OF  EFFECTIVE  THERMOVISCOELASTIC  PARAMETERS 


Note  that  in  the  equations  of  Section  2.1,  the  stresses  and 
temperatures  are  applied,  and  the  composite  or  average  strains  are  assumed 
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to  be  known  so  that  the  effective  compliance  parameters  may  be  calculated. 
The  analytical  tools  available  for  generating  average  strains  in  composite 
cylinders  involving  transversely  isotropic  fibers  and  an  isotropic  matrix 
include  the  composite  cylinder  assemblage  (CCA)  model  and  the  analysis  of  a 
periodic  hexagonal  array  [17].  The  details  regarding  the  analysis  for 
average  strains  in  a  composite  cylinder  composed  of  transversely  isotropic 
temperature -dependent  elastic  fibers  and  a  thermorheologically  complex 
viscoelastic  matrix  were  developed  in  the  Phase  I  portion  of  this  effort 
[16] .  Average  strains  due  to  thermal  expansion  and  axisymmetric  mechanical 
loadings  were  examined  with  the  composite  cylinder  assemblage  model,  while 
transverse  loading  effects  were  analyzed  using  a  finite  element  model  of  a 
periodic  hexagonal  array. 

Numerical  simulations  of  the  series  of  isothermal  creep  tests 
described  in  Section  2.1  were  performed  in  each  mode  of  deformation,  to 
determine  the  complete  set  of  effective  creep  compliance  parameters  of  the 
composite  constitutive  equations.  Figures  2.1  through  2.3  show  the  time 
dependent  creep  compliance,  horizontal  shift  and  vertical  shift, 
respectively,  for  the  composite  axial  deformation  compliance.  Since  this 
deformation  mode  is  largely  dominated  by  the  elastic  fibers,  the  creep 
compliance  (Figure  2.1)  is  several  orders  of  magnitude  less  than  that  of  the 
matrix  shear  creep  compliance.  This  fiber  dominance  also  manifests  itself 
in  the  shifts  (Figures  2.2  and  2.3)  through  a  much  more  shallow  slope  in  the 
composite  shift  curves  than  in  the  matrix  shift  curves.  Figures  2.4  through 
2.6  show  the  time  dependent  creep  compliance,  horizontal  shift  and  vertical 
shift,  respectively,  for  the  composite  axial  shear  compliance.  Since  in 
this  mode,  the  matrix  participation  is  dominant,  the  composite  parameters 
are  much  closer  to  the  matrix  values. 


2.4.  VERIFICATION  OF  EFFECTIVE  COMPOSITE  CONSTITUTIVE  EQUATIONS 

The  macromechanical  response  of  a  composite  structural  element,  as 
predicted  by  the  effective  constitutive  equations  and  their  derived 
parameters,  was  checked  against  results  computed  using  a  micromechanical 
model  (composite  cylinder  assemblage)  which  explicitly  included  the  fiber 
and  matrix  as  discrete  phases.  A  slowly  varying  periodic  shear  stress  was 
applied  to  the  composite  in  the  presence  of  an  oscillating  temperature.  It 
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should  be  noted  that  even  though  the  composite  cylinder  assemblage  model  was 
used  to  generate  the  parameters  of  the  effective  thermoviscoelastic 
equations,  this  model  can  still  be  used  to  verify  the  equations  and  their 
parameters  for  arbitrary  time  varying  mechanical  and  temperature  loadings, 
since  the  parameters  were  generated  using  a  series  of  isothermal  creep 
tests . 

The  simultaneously  applied  axial  stress  history  and  temperature  load 
are  shown  in  Figures  2.7  and  2.8,  respectively.  Comparisons  of  the  axial 
shear  strains  and  the  thermally  induced  longitudinal  and  transverse  strains 
for  this  loading  are  shown  in  Figures  2.9,  2.10  and  2.11  respectively.  Note 
that  essentially  identical  results  for  all  of  the  strains  were  obtained  from 
the  macro-  and  micromechanical  models,  thereby  verifying  both  the  assumed 
form  and  the  individual  parameters  of  the  composite  constitutive  equations. 

The  significance  of  the  verification,  illustrated  in  Figures  2.9 
through  2.11,  is  that  it  allows  structural  analyses  involving  TCM  composites 
to  be  performed  in  an  efficient  manner,  namely  by  simply  using  the  effective 
constitutive  equations.  If  these  equations  could  not  be  verified,  the 
response  of  the  composite  element  would  have  to  be  computed  at  each  time 
step  by  using  both  the  CCA  model  and  the  finite  element  model  of  the 
hexagonal  array,  both  of  which  require  the  solution  of  simultaneous 
algebraic  equations,  just  to  evaluate  the  average  strain  field  in  the 
composite  due  to  the  applied  mechanical  and  thermal  loads.  This  would  be  an 
extremely  tedious  and  computationally  excessive  process  for  any  practical 
structural  analysis  involving  TCM  composites. 
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3.  RESPONSE  TO  SINUSOIDAL  TEMPERATURE  AND  LOAD 


In  isothermal  viscoelasticity,  the  application  of  a  sinusoidal 
stress  results  in  a  sinusoidal  strain  at  the  same  frequency.  The 
oscillating  stress  and  strain  reach  their  peak  values  at  different  times, 
i.e.  there  is  a  phase  shift  between  stress  and  strain.  This  can  be  seen  in 
Figures  3.1  and  3.2  which  show  respectively  a  sinusoidal  axial  shear  stress 
history  applied  isothermally  and  the  resulting  axial  shear  strain.  The 
strain  response  is  most  conveniently  characterized  using  a  complex 
compliance,  where  the  real  part  is  the  ratio  of  the  strain  amplitude  to 
stress  amplitude,  and  the  imaginary  part  represents  the  phase  shift.  At  the 
frequency  of  the  applied  stress  shown  in  Figure  3.1,  both  the  real  and 
imaginary  parts  of  the  complex  compliance  can  be  computed  using  the 
information  on  the  peak  magnitudes  and  the  time  difference  between  the  peaks 
shown  in  Figures  3.1  and  3.2.  In  the  general  case,  both  the  real  and 
imaginary  parts  are  functions  of  the  frequency  of  the  applied  stress. 

The  purpose  of  analyzing  the  effects  of  simultaneous  sinusoidal 
temperature  and  load  was  to  determine  if  a  periodic  response  was  obtained 
for  the  composite,  and  if  so,  to  examine  the  form  of  the  periodic  response. 
In  particular  it  was  of  interest  to  determine  if  a  complex  compliance 
existed  in  this  situation. 

Several  cases  were  run  involving  the  application  of  a  sinusoidal 
shear  stress  and  a  sinusoidal  temperature  at  different  frequencies.  Figure 

3.3  shows  the  resulting  axial  shear  strain  for  the  case  where  the  shear 
stress  frequency  was  twice  the  frequency  of  the  temperature  cycle.  Figure 

3.4  gives  the  shear  strain  for  the  case  where  the  shear  stress  frequency  was 
two -thirds  the  applied  temperature  frequency,  and  Figure  3.5  gives  the  shear 
strain  when  the  shear  stress  frequency  is  two- fifths  the  applied  temperature 
frequency.  It  is  apparent  from  all  of  these  curves  that  the  steady  state 
portion  of  the  strain  can  be  decomposed  into  two  sinusoidal  components, 
whose  frequencies  are  equal  to  the  applied  stress  frequency  and  the 
temperature  frequency.  This  implies  the  existence  of  a  two  frequency 
complex  compliance  which,  from  a  computational  point  of  view,  could  be 
developed  into  a  powerful  tool  for  performing  structural  dynamic  analysis  of 
TSM  and  TCM  structures  under  arbitrary  mechanical  and  thermal  loading 
conditions . 


4.1.  INTRODUCTION 


The  purpose  of  this  section  is  to  present  the  results  of  the  effects 
of  thermoviscoelastic  behavior  on  the  damped  response  of  some  simple 
composite  structures.  Two  separate  problems  were  considered:  the  free 
vibration  of  a  cantilevered  unidirectional  beam  of  cylindrical  cross- 
section,  and  the  response  of  a  simple  truss  to  impulsive  mechanical  loading 
in  the  presence  of  isothermal  and  transient  temperature  conditions.  The 
role  of  shear  deformation  is  demonstrated  in  the  free  vibration  study,  and 
the  influence  of  matrix  deformation  on  damped  response  is  shown  in  the  truss 
problems . 

4.2.  FREE  FLEXURAL  VIBRATIONS  OF  CANTILEVERED  UNIDIRECTIONAL  BEAM 

As  an  initial  illustration  of  dynamic  analysis  of  a  TCM  composite 
structure,  we  consider  the  case  of  free  flexural  vibrations  of  a 
«  cantilevered  beam  which  is  uniaxially  reinforced  in  the  beam  axis  direction. 

As  mentioned  above,  the  purpose  of  the  investigation  is  to  compare  vibration 
damping  due  to  matrix  thermoviscoelasticity  on  the  basis  of  the  usual  theory 
which  neglects  the  effect  of  shear,  and  on  the  basis  of  the  more  refined 
Timoshenko  theory  which  takes  into  account  shear  as  well  as  rotatory 
inertia.  Recall  that  the  loss  tangent  is  defined  as  the  ratio  of  the 
imaginary  part  of  the  complex  modulus  to  its  real  part.  For  isotropic 
materials,  in  which  the  complex  Young's  modulus  loss  tangent  and  the  complex 
shear  modulus  loss  tangent  are  of  the  same  order,  the  added  effect  of  shear 
and  rotatory  inertia  is  small  for  vibration  modes  of  low  order  and  for  long 
beams.  In  the  present  case,  however,  where  the  axial  Young's  modulus  loss 
tangent  is  by  an  order  of  magnitude  smaller  than  that  of  the  axial  shear 
modulus,  the  situation  is  quite  different  as  will  be  shown. 

Considering  only  the  effect  of  flexure,  by  the  correspondence 
principle  for  viscoelastic  vibrations  [18],  the  differential  equation  for 
the  transverse  deflection  y  of  a  freely  vibrating  beam  is 
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(4.2.1) 
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pA  ^  - 
dz 


where 

—  R  I 

E  (u>)  -  E  (w)  +  iE  (w)  -  complex  axial  Young's  modulus, 

AAA 

I  -  moment  of  inertia 
A  —  area  of  the  cross  section 
p  -  density. 

By  the  conventional  separation  of  variables  method,  the  solution  to 
the  governing  partial  differential  equation  (4.2.1)  is  given  by 

y(x,t)  -  (A^  sin  ax  +  A2COS  ax  +  A^  sinh  ax  +  A^  cosh  ax)elut  (4.2.2) 

with 


4 

a 


pAD2/E  I 

A 


(4.2.3) 


where  the  constants  A^  are  evaluated  from  boundary  conditions  on  a 
cantilevered  beam 


y(0,t)  -  0 

y'(0,t)  -  0 

(4.2.4) 


M(L.t)  -  EAIy"(L,t)-0 
V(L,t)  -  EAIy"'(L,t)  -  0 
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Substituting  the  solution  given  by  (4.2.2)  into  the  homogeneous  boundary 
conditions,  by  requiring  a  non-trivial  solution  for  at  least  one  of  the 
constants,  leads  to  the  frequency  equation 


1  +  cos  aL  cosh  aL  -  0 


(4.2.5) 


The  solution  to  this  transcendental  equation  will  be  represented  as  a^L  so 


that  the  beam  natural  frequencies  may  be  obtained  from  the  solution  of 


Z  -  a2»(EA(w)I/pA)1/2 
n  n  A 


(4.2.6) 


Note  that  in  the  elastic  case,  the  natural  frequencies  are  obtained  by 
simply  evaluating  the  right  hand  side  of  (4.2.6).  However  in  the 


viscoelastic  case,  since  EA  is  a  function  of  frequency,  equation  (4.2.6)  is 


nonlinear  in  frequency  and  thus  the  real  and  imaginary  parts  must  be 
obtained  in  an  iterative  fashion. 

The  time  varying  portion  of  the  solution  takes  the  form 
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(4.2.7) 


Thus  the  real  part  of  each  natural  frequency  represents  the  oscillatory 
component  and  the  imaginary  part  is  the  attenuation,  a  measure  of  the 
damping  introduced  by  the  viscoelastic  matrix  material. 

Next  the  same  beam  is  considered  in  Timoshenko  fashion,  with  shear 
and  rotatory  inertia.  Again  using  the  correspondence  principle,  the  coupled 
equations  for  the  total  deflection  y  and  bending  slope  \J>  are  given  by 
Timoshenko  [19]  as 
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where  Che  new  constants  are 


—  R  I 

Ga(«)  ”  GA(w)  +  iG^(w)  ”  the  axial  complex  shear  modulus,  and 


k  -  a  numerical  shape  factor  for  the  cross  section. 


Following  the  development  first  shown  in  Huang  [20] ,  we  eliminate  ^  or  y 
from  equations  (4.2.8)  and  (4.2.9)  and  obtain  the  two  uncoupled  differential 
equations  in  y  and  ip 


au  y 


a1  y 


eax  "~4  +  pk  — 9  '  pI 

dx  at ^ 


1  + 


kG. 


4  4 

a  y  2T  a  y 

— -  +  e-L  —L  _  0 

ax2at2  kG  at4 

A 


(4.2.10) 
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Introducing  the  nondimens ional  parameter  9  -  x/L  and  assuming  solutions  to 
be  of  the  form 


y  -  Ye 
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(4.2.12) 


equations  (4.2.8)  to  (4.2.11)  take  the  form 


S2*"  -  (l-B2r2S2)*  +  Y'/L  -  0 


(4.2.13) 


(4.2.14) 


K  where 

8  c2  .4-2 

s  b  '  V  L " 

(4.2.17) 

y  2  2 

;5  r  -  1/ Al 

(4.2.18) 

|  S2  -  EAI/kAGAL2 

(4.2.19) 

o 

i 

(4.2.20) 

o 

I 

>* 

»N 

and  at  the  free  end  we  have 
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where 


—  Cc  cosh  ba#  +  C,  sinh  ba#  +  C-,  cos  b B8  +  CQ  sin  b 8S  (4.2.23 

J  O  /  O 


a,  0 


-^{^(^[(AiV^/tV'2}  1/2 


(4.2.24) 


Note  that  only  half  of  the  constants  in  equations  (4.2.22)  and 
(4.2.23)  are  independent.  The  constants  C,.  through  Cg  may  be  related  to 

through  using  equations  (4.2.13)  or  (4.2.14). 

From  the  requirement  that  a  non- trivial  solution  exist  for  at  least 
one  of  the  constants,  using  the  appropriate  boundary  conditions  (4.2.20)  and 
(4.2.21)  within  (4.2.22)  and  (4.2.23)  we  get  as  a  frequency  equation 


2+[b2(r2-S2)2+2]cosh  ba  cos  bp  -  — 3-C.F-+S2) 


(l-S2r2S2)1/2 


sinh  ba  sin  ba  -  0  (4.2.25) 


As  with  (4.2.6)  above,  a  nonlinear  solution  technique  must  be  used  on 
(4.2.25)  to  obtain  the  real  and  imaginary  components  of  the  beam  natural 
frequencies . 

To  use  the  natural  frequency  equations  given  above,  the  effective 
complex  Young's  modulus  and  the  effective  complex  axial  shear  modulus  for 
the  unidirectional  composite  must  be  known.  These  were  generated  in  a 
manner  entirely  analogous  to  that  used  in  Section  2  in  generating  the 
effective  creep  compliance  parameters  of  the  composite.  In  this  case, 
however,  instead  of  dealing  with  "reduced  time"  as  the  independent  variable 
of  the  time  dependent  creep  compliances,  we  work  with  "reduced  frequency"  as 
the  abscissa  of  the  complex  moduli.  As  with  the  creep  compliances,  the 
horizontal  shifts  are  used  to  compute  the  reduced  frequency  for  the  complex 
moduli . 

As  an  illustration  of  the  dependence  of  response  on  temperature,  the 
natural  frequencies  of  a  fixed- free  carbon-epoxy  unidirectional  tube  40 
inches  long,  0.01  inch  in  thickness,  and  1  inch  in  diameter  were  obtained  by 


solving  the  frequency  equations  (4.2.6)  and  (4.2.25)  under  isothermal 
conditions.  As  with  the  earlier  examples,  T300  fibers  were  embedded  with  a 
volume  fraction  of  0.6  within  a  Hercules  3502  epoxy  resin  matrix. 

Figure  4.1  shows  the  first  mode  natural  frequencies  for  the  carbon- 
epoxy  cantilevered  beam  as  a  function  of  temperature.  This  curve  shows  that 
the  real  part,  or  resonant  frequency,  is  essentially  independent  of 
temperature  and,  in  addition,  is  practically  unchanged  by  the  presence  of 
shear  deformation  and  rotatory  inertia.  The  imaginary  or  attenuation 
component  varies  considerably  with  temperature,  however,  and  is  strongly 
affected  by  shear  deformation.  The  effects  of  shear  deformation  can  be 
readily  seen  in  Figure  4.2,  which  shows  the  envelopes  of  the  first  mode 
response  of  both  beam  models  at  room  temperature.  Figure  4.3  shows  the 
second  mode  natural  frequencies  for  the  beam  as  a  function  of  temperature. 
For  the  Bernoulli-Euler  beam  the  resonant  frequency  remains  independent  of 
temperature,  but  both  shear  deformation  and  temperature  have  an  effect  on 
the  second  resonant  frequency,  as  can  be  seen  by  the  difference  occurring  in 
the  real  components.  Shear  deformation  effects  on  the  imaginary  part  are 
even  more  pronounced  than  in  the  first  mode,  as  can  be  seen  by  the  increased 
spread  between  the  two  imaginary  components.  The  result  of  shear 
deformation  on  the  second  mode  is  even  more  apparent  in  Figure  4.4,  which 
shows  a  significantly  shorter  decay  time  for  the  response  of  the  Timoshenko 
beam  than  for  the  Bernoulli-Euler  beam. 

4.3.  EFFECT  OF  TRANSIENT  LOADS  AND  TEMPERATURES  ON  A  SIMPLE  COMPOSITE 

STRUCTURE 

As  a  second  illustration  of  dynamic  analysis  of  a  TCM  composite 
structure,  we  consider  the  case  of  the  simple  truss  shown  in  Figure  4.5. 

All  three  truss  members  are  unidirectional  carbon-epoxy  composite  rods, 
characterized  using  TCM  behavior.  Note  from  this  figure  that  a  zero  degree 
layup  is  defined  as  fibers  parallel  to  the  axis  of  the  bar  and  a  90  degree 
layup  is  defined  as  fibers  perpendicular  to  the  axis  of  the  bar. 

To  perform  the  truss  dynamic  analyses,  rod  finite  elements  were 
developed  using  the  standard  linear  displacement  field  [23]  for  uniaxial 
members  composed  of  thermorheologically  complex  material.  The  TCM  rod 
elements  were  then  incorporated  within  a  finite  element  code  capable  of 


performing  isothermal  natural  frequency  analysis,  as  well  as  static, 
thermal,  and  dynamic  stress  analysis.  The  dynamic  analyses  were  performed 
in  the  time  domain  using  the  direct  integration  procedure  of  Newmark  [24]  as 
outlined  by  Bathe  and  Wilson  [25]. 

As  a  first  step  towards  examining  truss  response  under  transient 
loading  conditions,  an  eigensolution  of  the  two  degree-of-freedom  truss  was 
performed  at  a  temperature  of  117  degrees  Kelvin  above  the  reference 
temperature  (303  degrees  Kelvin) .  The  mode  shapes  and  natural  frequencies 
resulting  from  the  eigensolution  are  given  in  Figure  4.6.  The  natural 
frequencies  for  a  truss  composed  entirely  of  zero  degree  layup  members  and  a 
truss  composed  entirely  of  90  degree  layup  members  are  shown  in  the  figure. 

The  truss  was  subjected  to  an  impulsive  horizontal  load  1  kN  in 
magnitude  and  one  second  in  duration,  and  the  response  was  examined  for 
isothermal  as  well  as  transient  temperature  conditions  applied  uniformly 
over  the  entire  truss.  At  117  degrees  Kelvin  above  the  reference 
temperature,  Figure  4.7  shows  the  vertical  displacement  at  the  apex  (joint 
2)  of  the  truss  with  zero  degree  layup  members.  In  this  truss,  most  of  the 
load  is  carried  by  the  fibers,  so  that  there  is  little  matrix  deformation 
occurring  during  the  response.  This  results  in  very  little  damping  as  shown 
by  the  insignificant  amount  of  decay  occurring  in  the  one  minute  of  response 
shown  in  the  figure.  For  the  same  loading  and  at  the  same  temperature,  the 
vertical  displacement  at  joint  2  of  the  truss  with  90  degree  layup  members 
is  shown  in  Figure  4.8.  In  addition  to  the  larger  response  period, 
resulting  from  90  degree  layup  members  having  smaller  stiffness  than  zero 
degree  layup  members,  a  substantial  amount  of  decay  in  the  response  is 
evident  from  this  figure,  a  direct  result  of  more  matrix  deformation 
occurring  in  this  configuration  than  in  the  zero  degree  layup  configuration. 

For  the  same  impulsive  mechanical  load,  two  other  thermal  load 
conditions  were  examined  for  the  truss  with  90  degree  layup  members.  The 
response  at  the  reference  temperature  as  well  as  the  response  occurring  for 
the  transient  temperature  condition  of  A#  -  117°K  x  t  (min)  were  both 
examined.  For  these  two  load  cases,  peak  responses  were  between  the  limits 
of  the  peaks  of  Figures  4.7  and  4.8.  The  envelopes  of  the  vertical 
displacement  response  histories  for  all  four  cases  examined  are  shown  in 
Figure  4.9.  Not  surprisingly,  as  in  the  cantilever  free  vibration  problem, 


the  amount  of  damping  present  in  the  response  is  dependent  on  temperature 
and  the  extent  of  matrix  deformation  in  the  composite  members. 

As  a  final  truss  analysis  problem,  the  response  of  the  truss 
composed  entirely  of  zero  degree  layup  members  and  with  the  lumped  mass  at 
joint  2  removed,  was  examined  under  thermal  shock  conditions,  in  the  absence 
of  any  mechanical  loading.  A  temperature  varying  linearly  from  80°K  to 
420°K  in  0.5  seconds  and  held  constant  at  420°K  thereafter  was  used  as  the 
thermal  load  applied  uniformly  to  all  members  of  the  truss,  as  shown  in 
Figure  4.10.  The  horizontal  and  vertical  displacement  histories  at  joint  2 
of  the  truss  are  shown  in  Figures  4.11  and  4.12  respectively.  Due  to 
changes  in  the  displacement  and  velocity  occurring  rapidly  enough  for 
inertia  forces  to  be  significant,  oscillations  occur  during  the  entire 
response  interval  examined  and,  when  the  temperature  becomes  constant  at  0.5 
seconds,  the  oscillations  take  place  about  the  displacement  reached  at  that 
^ime.  The  stresses  in  members  1  and  2  are  shown  in  Figures  4.13  and  4.14 
respectively.  Since  joint  2  is  unrestrained,  these  stresses  result  only 
from  the  oscillations  in  the  displacements  occurring  at  joint  2. 


DISCUSSION  AND  RECOMMENDATIONS 


The  research  program  detailed  in  this  report  has  had  several 
significant  results,  among  them  the  derivation  and  verification  of  the 
composite  TCM  stress-strain  equations,  and  the  demonstration  of  the 
significance  of  temperature  and  shear  deformation  on  the  damped  response  of 
TCM  composite  structures.  From  the  standpoint  of  performing  dynamic 
analysis  on  space  structures  composed  of  carbon  fiber/polymeric  matrix 
composites,  the  results  of  applying  sinusoidal  mechanical  and  thermal  loads 
on  TCM  composite  structural  elements,  and  observations  concerning  the 
computations  in  the  forced  vibration  analysis  of  the  truss  structures  were 
equally  important. 

As  discussed  in  Section  3,  numerical  experiments  performed  indicate 
the  existence  of  two  frequency  complex  compliances/moduli  for 
thermorheologically  complex  materials .  By  analogy  to  isothermal 
viscoelastic  analysis,  this  means  that  for  periodic  thermal  and  mechanical 
loads,  the  TCM  response  will  be  most  conveniently  obtained  by  the 
superposition  of  the  effects  due  to  the  applied  stress  frequency  components 
plus  the  thermal  load  frequency  components.  In  other  words,  for  periodic 
loads,  frequency  domain  response  is  more  efficient. 

The  dynamic  analyses  of  the  truss  structure  detailed  in  Section  4 
were  accomplished  by  direct  time  integration  of  the  differential  equations 
of  motion.  At  each  time  step  and  for  each  element,  the  hereditary  integrals 
must  be  re-evaluated.  In  the  case  of  TCM  creep  compliances  represented 
using  power  law  forms,  this  is  a  computationally  intensive  task  which  takes 
more  and  more  effort  as  the  duration  of  response  becomes  larger  and  larger. 
Consequently  for  isothermal  truss  dynamic  analyses  involving  either  periodic 
or  non-periodic  loads,  the  complex  frequency  response  method  is  far  more 
computationally  efficient.  With  the  development  of  the  two  frequency 
complex  compliances/moduli  for  TCM  members,  the  complex  frequency  response 
method  could  be  extended  to  allow  the  solution  of  periodic  and  non-periodic 
thermal  and  mechanical  loads  to  be  performed  in  the  frequency  domain.  In 
addition,  the  versatility  gained  by  using  a  frequency  domain  formulation,  in 
particular  the  ability  to  perform  stochastic  analyses,  would  enhance  the 
attractiveness  of  this  approach. 
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Table  1.  Material  Properties 


Fiber  Properties 

Ea  -  231  GPA 
A 

Ex  -  22.4  GPa 

Ga  =  22.1  GPa 
A 

Gx  -  8.30  GPa 

u.  =  0.30 
A 

a =  -1.33.  pm/m°c 

=  7.04  /iin/m°c 

VD  =  0.60 
F 

Matrix  properties  at  reference  temperature  (303°K) 

k  -  3.58  GPa 

g(o)  -  0.606  (GPa)'1 

i/  -  0.30 

a  =>30.0 


35 


&&&&& 14 


<  i/Pa) 


1  2  3 

Time  (  hours  ) 


Curve  Legend 


Matrix  Shear  Compliance 
Composite  11  Compliance 


Figure  2.1.  Effective  Axial  Deformation  Creep  Compliance  for  Unidirectional 
Carbon/Epoxy  Composite 
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Figure  2.3.  Effective  Axial  Deformation  Vertical  Shift  for  Unidirectional 
Carbon/Epoxy  Composite 
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Effective  Axial  Shear  Creep  Compliance  for  Unidirectional 
Carbon/Epoxy  Composite 
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Figure  2.5,  Effective  Axial  Shear  Horizontal  Shift  for  Unidirectional 
Carbon/Epoxy  Composite 


WWW 


TiMe  (  hours  ) 


*10* 


Figure  2.7.  Applied  Axial  Shear  Stress  History  on  Unidirectional 
Carbon/Epoxy  Composite  Structural  Element 
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Figure  2.9.  Comparison  of  Axial  Shear  Strain  Results  for  Simultaneous 
Axial  Shear  Stress  and  Temperature  Loads 
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Figure  2.10.  Comparison  of  Thermally  Induced  Longitudinal  Strains 
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Figure  2.11.  Comparison  of  Thermally  Induced  Transverse  Strains 
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Figure  4.4.  Envelope  of  Second  Mode  Response  in  Unidirectional  Composite  Cantilever 
Beam  at  Room  Temperature 


L  *»  20  m,  M  =  870  kg,  Carbon-Epoxy  Truss  Members 


0°  Layup  -  Fibers  Parallel  to  Axis  of  Bar 

90°  Layup  -  Fibers  Perpendicular  to  Axis  of  Bar 


Figure  4.5.  Geometry  of  Simple  Truss  Used  in  Composite  Structural 
Analyses 
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Figure  4.6.  Eigensolution  of  Simple  Truss  with  Zero  and  90 
Degree  Layup  Composite  Members 
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Figure  4.7.  Vertical  Displacement  at  Joint  2  Due  to  Impulsive 
Horizontal  Load  at  A4>  =  117°K  for  Truss  with  Zero 
Degree  Layup  Members 
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Figure  4.8.  Vertical  Displacement  at  Joint  2  Due  to  Impulsive 
Horizontal  Load  at  At  =  117°K  for  Truss  with  90 
Degree  Layup  Members 
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.12.  Vertical  Displacement  History  at  Joint  2  Due  to 
Thermal  Shock  Load 
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Figure  4.13.  Member  1  Stress  History  for  Applied  Thermal  Shock 
Load 


